L-selectin mediated leukocyte tethering in shear flow is controled by multiple contacts 
and cytoskeletal anchorage facilitating fast rebinding events 
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L-selectin mediated tethers result in leukocyte rolling only above a threshold in shear. Here we 
present biophysical modeling based on recently published data from flow chamber experiments (Dwir 
et al., J. Cell Biol. f63: 649-659, 2003) which supports the interpretation that L-selectin mediated 
tethers below the shear threshold correspond to single L-selectin carbohydrate bonds dissociating 
on the time scale of milliseconds, whereas L-selectin mediated tethers above the shear threshold are 
stabilized by multiple bonds and fast rebinding of broken bonds, resulting in tether lifetimes on the 
timescale of 10 _1 seconds. Our calculations for cluster dissociation suggest that the single molecule 
rebinding rate is of the order of 10 4 Hz. A similar estimate results if increased tether dissociation 
for tail-truncated L-selectin mutants above the shear threshold is modeled as diffusive escape of 
single receptors from the rebinding region due to increased mobility. Using computer simulations, 
we show that our model yields first order dissociation kinetics and exponential dependence of tether 
dissociation rates on shear stress. Our results suggest that multiple contacts, cytoskeletal anchorage 
of L-selectin and local rebinding of ligand play important roles in L-selectin tether stabilization and 
progression of tethers into persistent rolling on endothelial surfaces. 



Introduction 



Leukocyte traffiking plays a central role in the immune 
response of vertebrates. Leukocytes constantly circulate 
in the cardiovascular system and enter into tissue and 
lymph through a multi-step process involving rolling on 
the endothelium, activation by chemokines, arrest, and 
transendothelial migration A key molecule in this 
process is L-selectin, a leukocyte-expressed adhesion re- 
ceptor which is localized to tips of microvilli and binds 
to glycosylated ligands on the endothelium. Its proper- 
ties are optimized for initial capture and rolling under 
physiological shear 0j , as confirmed byrecent exper- 
imental data and computer simulations |3, In con- 
trast to tethering through other receptor systems like 
P-selectin, E-selectin or integrins, appreciable tether- 
ing through L-selectin and subsequent rolling only oc- 
curs above a threshold in shear [(|, even in cell- free sys- 
tems 0, • Downregulation by low shear is unique for 
L-selectin tethers and might be necessary because L- 
selectin ligands are constitutively expressed on circulat- 
ing leukocytes, platelets and subsets of blood vessels 0- 

The dissociation rate of single molecular bonds is ex- 
pected to depend exponentially on an externally applied 
steady force (Bell equation) |10| . Quantitative analy- 
sis with regular video camera (time resolution of 30 ms) 
of L-selectin tether kinetics in flow chambers above the 
shear threshold resulted in first-order dissociation kinet- 
ics, with a force dependence which could be fit well to the 
Bell equation, resulting in a force-free dissociation con- 
stant of 6.6 Hz 0, H |1 H3. These findings have been 
interpreted as signatures of single L-selectin carbohy- 
drate bonds. However, recent experimental evidence sug- 
gests that L-selectin tether stabilization involves multiple 



bonds and local rebinding events. Evans and coworkers 
used the biomcmbrane force probe to measure unbinding 
rates for single L-selectin bonds as a function of loading 
force ^2 • Modeling bond rupture as thermally activated 
escape over a sequence of transition state barriers in- 
creasingly lowered by rising force [U| , these experiments 
revealed two energy barriers along the unbinding path- 
way. The inner barrier corresponds to Ca 2+ -dependent 
binding through the lectin domain and explains the high 
strength of L-sclcctin mediated tethers required for cell 
capture from shear flow. Extracting barrier properties 
from the dynamic force spectroscopy data allows to con- 
vert them into a plot of dissociation rate as a function 
of force. In this way, results from dynamic force spec- 
troscopy and flow chamber experiments can be compared 
in a way which is independent of loading rate. In detail, 
Evans and coworkers found a 1000-fold increase in disso- 
ciation rate as force rises from to 200 pN, in marked 
contrast to tether dissociation kinetics as measured in 
flow chamber experiments, which increases at most 10- 
fold over this range 0, Q ■ Therefore additional stabiliza- 
tion has to be involved with leukocyte tethers mediated 
by L-selectin. Dwir and coworkers used flow chambers 
to study tethering of leukocytes transfected with tail- 
modified mutants of L-selectin [14j|. They found that 
tether dissociation increases with increased tail trunca- 
tion, possibly since tail truncation leads to decreased cy- 
toskeletal anchorage and increased mobility. More re- 
cently, Dwir and coworkers found with high speed video 
camera (time resolution of 2 ms) that L-selectin teth- 
ers form even below the shear threshold at shear rate 
40 Hz, albeit with a very fast dissociation rate of 250 
Hz, undetectable with regular camera ^5|. Thus the 
shear threshold results from insufficient tether stabiliza- 
tion at low shear. Using systematic changes in viscosity 
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(which changes shear stress, but not shear rate), Dwir 
and coworkers were able to show that at the shear thresh- 
old, tether lifetime is prolongued by a factor of 14 due 
to shear-mediated cell transport over L-selectin ligand. 
They suggested that sufficient transport might be needed 
for formation of additional bonds. With more than one 
bond being present, rebinding then could provide the 
tether stabilization observed experimentally. 

In this paper, we present a theoretical model for the 
interplay between bond rupture, L-selectin mobility and 
ligand rebinding within small clusters of L-selectin bonds, 
which interprets the recent experimental results in a con- 
sistent and quantitative way. Traditionally, tether dis- 
sociation at low ligand density has been interpreted as 
single molecule rupture due to the observed first order 
dissociation kinetics and a shear dependence which can 
be fit well to the Bell equation. Here we demonstrate 
that the same features result for small clusters of multi- 
ple bonds with fast rebinding. Our results suggest that 
the shear threshold corresponds to the formation of mul- 
tiple contacts and that single L-selectin bonds decay too 
rapidly as to provide functional leukocyte tethers. 



Experiments 

Our experimental procedures have been described be- 
fore elsewhere [hH Il5j . Three variants of human L- 
selectin were stably expressed in the mouse 300.19 pre 
B cell line. Wildtype, tail-truncated mutant and tail- 
deleted mutant have the same extracellular domains and 
differ only in their cytoplasmic tails. L-selectin medi- 
ated tethering was investigated in a parallel plate flow 
chamber. The main ligand used was PNAd, the major 
L-selectin glycoprotein ligand expressed on endothelium. 
For immobilization in the flow chamber, the ligand was 
diluted in such a way that no rolling was supported at 
shear rates lower than 100 Hz (dilution 10 ng/ml in the 
coating solution, which corresponds to an approximate 
scaffold density of 100//im 2 ). Single tethers were moni- 
tored with video microscopy at 2 ms resolution and the 
microkinetics were analyzed by single cell tracking as de- 
scribed previously. The logarithm of the number of cells 
which pause longer than time t is plotted as a function 
of t and usually gives a straight line indicative of an ef- 
fectively first order dissociation process. The slope is the 
tether dissociation rate k ff and is plotted as a func- 
tion of shear rate 7 in Fig. 1. This plot shows that be- 
low the shear threshold of 40 Hz, the dissociation rate 
is 250 Hz, independent of tail mutations and viscosity 
of the medium. Above the shear threshold, the dissoci- 
ation rate becomes force-dependent, with a dependence 
on shear stress which can be fit well to the Bell equation 
k a ff = fc e F / Fb 10]. Here fc is the force- free dissociation 
rate and Fj, is the bond's internal force scale. The force 
on an undeformed 300.19 lymphocyte with radius R = 6 
fim. follows from Stokes flow around a sphere close to 
a wall [16|. Taking into account the lever arm geometry 




(b) (c) 

FIG. 1: Tether dissociation rate k ff determined from kinetic 
analysis of flow chamber experiments plotted as function of 
shear rate 7 |T^. Solid line with circles: wildtype. Dashed 
line with diamonds: tail-deleted mutant. Dotted line with 
squares: wildtype with 6 percent of Ficoll, which changes 
viscosity and thus shear stress (but not shear rate) by a fac- 
tor of 2.6. These data suggest that the shear threshold is a 
transport-related rather than a force-related issue, and that 
the shear threshold is not about ligand recognition, but about 
tether stabilization. 



provided by the tether holding the cell at an angle of 50°, 
the force acting on the L-selectin bonds can be calculated 
to be F = 180 pN per dyn/cm 2 of shear stress Fitting 
the Bell equation to the wildtype data from Fig^ 1 give s 
similar values as obtained in earlier studies 
namely a force- free dissociation constant of 6.6 Hz and an 
internal force scale Fj, = 200 pN (corresponding to a reac- 
tive compliance of 0.2 A). At the shear threshold, we find 
14-fold and 7-fold reduction in dissociation rate for wild- 
type and tail-truncated mutant, respectively. Adding 6 
volume percent of the non-toxic sugar Ficoll increases vis- 
cosity from 1 cP to 2.6 cP. Thus shear stress is increased 
by a factor of 2.6, whereas shear rate is unchanged. At 
the shear threshold, this increases wildtype dissociation 
3-fold, roughly as expected from the fit to the Bell equa- 
tion. Most importantly, there is no shift of the shear 
threshold as a function of shear rate. This indicates that 
the shear threshold results from shear-mediated trans- 
port, rather than from a force-dependent process. 
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FIG. 2: A schematic representation of the mechanisms in- 
volved in L-selectin mediated leukocyte tethering to diluted 
carbohydrate ligands. (a) Initial binding most likely corre- 
sponds to one L-selectin receptor localized to the tip of one 
microvillus binding to ligand presented on a glycoprotein scaf- 
fold on the substrate. At low shear, stabilization through 
additional bonds is unlikely, because the distance between 
scaffolds is larger than the microvilli's tips and the proba- 
bility of two microvilli simultaneously hitting two ligands is 
very low. (b) At sufficiently high shear, shear-mediated rota- 
tion of the cell over the substrate leads to the establishment of 
an additional bond on another microvillus. In contrast to this 
two-dimensional cartoon, in practice the two microvilli are ex- 
pected to coexist with similar latitude, so that they can share 
force in a cooperative way. (c) Close-up to the cell-substrate 
interface. The L-selectin receptor can move laterally in the 
membrane, with an effective diffusion constant which depends 
on cytoskeletal anchorage. If a receptor has bound to lig- 
and on the substrate, it will rupture in a stochastic manner, 
depending on shear-induced loading. If an additional bond 
(most likely on the second microvillus) holds the cell during 
times of rupture, rebinding can occur at the first microvillus, 
thus increasing tether stabilization. 



Theory 

Shear-mediated transport. At the shear thresh- 
old at shear rate 7 = 40 Hz (corresponding to shear 
stress t = 777 = 0.4 dyn/cm 2 for viscosity 77 = 1 cP) 
and for small distance between cell and substrate, a cell 
with radius R — 6 /im will translate with hydrodynamic 
velocity u = 0.48i?7 = 0.12 /jm/ms and at the same 
time rotate with frequency £1 = O.267 = 10.4 Hz [ilij . 
Therefore the cell surface and the substrate surface will 
move relatively to each other with an effective velocity 
v = u — RQ, = 0.22i?7 = 50 nm/ms. In average there is 
no normal force which pushes the cell onto the substrate, 
but since it moves in close vicinity to the substrate, it can 
explore it with this effective velocity v. Thus there ex- 
ists a finite probability for a chance encounter between a 



L-selectin receptors on the tip of a microvillus and a car- 
bohydrate ligand on the substrate. Here we focus on the 
case of diluted ligands, with a ligand density of 100//zm 2 . 
Then the average distance between single ligands is 100 
nm, that is larger than the lateral extension of the mi- 
crovilli, which is 80 nm. Therefore the first tethering 
event is very likely to be a single molecular bond (com- 
pare Fig. 2a). If this first bond has formed, the microvil- 
lus will be pulled straight and the cell will slow down. It 
will come to a stop on the distance x of order /jm (e.g. 
the rest length of a microvillus is 0.35 /im). This takes 
the typical time t s — x/u — 8 ms. During this time, 
the cell can explore an additional distance of the order 
of vt s = 400 nm. The experimental data presented in 
Fig. 1 suggest that this is the minimal transport required 
to establish a second microvillar contact which is able to 
contribute to tether stabilization (compare Fig. 2b). 

Single bond loading. If tether duration was much 
longer than the time over which the cell comes to a stop, 
the single bond dissociation rate k ff below the shear 
threshold should increase exponentially with shear rate 7 
according to the Bell equation. However, this assumption 
is not valid in our case, because tether duration and slow- 
ing down time are both in the millisecond range. Fig. 3 
shows that indeed the Bell equation (dotted line) does 
not describe the wildtype data from Fig. 1 (dashed line 
with circles). In order to model a realistic loading pro- 
tocol, we assume that the force on the bond rises linear 
until time t s and then plateaus at the constant force F 
arising from shear flow. Note that initial loading rate 
r = F/t s scales quadratically with shear rate 7, because 
F ~ 7 and t B ~ 1 /j- The dissociation rate k a f / for this 
situation can be calculated exactly. The result is given in 
the supplemental material and is plotted as dash-dotted 
line in Fig. 3. It is considerably reduced towards the ex- 
perimentally observed plateau. Agreement is expected to 
increase further if initial loading is assumed to be sub- 
linear. A scaling argument shows the main mechanism at 
work. For the case of pure linear loading, the mean time 
to rupture is T = (Fb/r) exp(koFjJr)E(koFb/r), where 
E(x) is the exponential integral ^21 • There are two dif- 
ferent scaling regimes for slow and fast loading, which 
are separated by the critical loading rate r c = fco-Ffc- For 
slow loading, r < r c , a large argument expansion gives 
T sw 1/ko, that is the bond decays by itself before it 
starts to feel the effect of force. For fast loading, a small 
argument expansion gives T w (Fb/r) hx(r/koFf,), which 
is also found for the most frequent time of rupture in 
this regime [ljj. In our case, k = 250 Hz, F b = 200 pN 
and r c — fco-Ffc = 5 x 10 4 pN/s. At the shear threshold, 
r = 10 4 pN/s and we are still in the regime of slow load- 
ing, r < r c . This suggests that tethers below the shear 
threshold correspond to single L-selectin carbohydrate 
bonds which decay before the effect of force becomes ap- 
preciable. This does not imply that the bonds do not feel 
any force (after all the cell is slowed down), but that we 
are in a regime in which k a f f as a function of shear does 
not change appreciably, as observed experimentally. 
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FIG. 3: Theoretical predictions for tether dissociation rate 
k ff as a function of shear rate 7 compared to experimentally 
measured wildtype data from Fig. 1 (dashed line with circles). 
Dotted line: the single bond dissociation rate with force-free 
dissociation rate fco = 250 Hz and constant instantaneous 
loading increases exponentially according to the Bell equation. 
Dash-dotted line: it is reduced towards the experimentally 
observed plateau below the shear threshold at 7 = 40 Hz 
by including the effect of finite loading rates. Solid lines from 
top to bottom: cluster dissociation rate for two-bonded tether 
with rebinding rate k on = 0, 10, 20,40 and 60 fco. Above the 
shear threshold, the two-bonded tether with k on « 10 4 Hz 
agrees well with the experimentally measured data. 



Single bond rebinding. Single bond rupture is a 
stochastic process according to the dissociation rate k f / 
given by the Bell equation. If ligand and receptor re- 
main in spatial proximity after rupture, rebinding be- 
comes possible. We define the single molecule rebinding 
rate k on to be the rate for bond formation when recep- 
tor and ligand are in close proximity. If bond formation 
was decomposed into transport-determined formation of 
an encounter complex and chemical reaction of the two 
partners, then k on would correspond to the on-rate for 
reaction jlOl Il8j . It has the dimension of 1 /s and should 
not be confused with two- or three-dimensional associa- 
tion rates, which have dimensions of m 2 /s (equivalently 
m/Ms) and m 3 /s (equivalently 1/Ms), respectively. k on 
should depend mainly on the extracellular side of the 
receptor. In the following, it will therefore be assumed 
to be the same for wildtype and mutants. There are 
two mechanisms which might prevent rebinding within 
an initially formed cluster: the single receptor might es- 
cape from the rebinding region due to lateral mobility, 
or the receptor might be carried away from the ligand 
because the cell is carried away by shear flow. Fig. 2c 
shows schematically the interplay between rupture, re- 
binding and mobility for single L-selectin receptors. 

We start with the first case, that is diminished re- 



binding due to lateral receptor mobility. Since increased 
tail truncation decreases interaction with the cytoskele- 
ton 01 > lateral mobility increases from wildtype through 
tail-truncated to tail-deleted mutant. For each receptor 
type, we assume an effective diffusion constant D. The 
conditional probability for rebinding depends on absolute 
time since rupture. We approximate it by the probabil- 
ity that a particle with two-dimensional diffusion, but 
without capture is still within a disc with capture ra- 
dius s at time t, k on (t) = k on (l — e~ s / 4Dt ). Thus the 
time scale for the diffusion correction is set by s 2 /AD, the 
time to diffuse the distance of the capture radius. The 
diffusion constant for the wildtype can be estimated to 
be 10 -11 cm 2 /s, with the one for the tail-deleted mutant 
being at least one order of magnitude larger [23, 0OI ■ A 
typical value for the capture radius is s = 1 nm. Then 
the time t c — s 2 /AD to diffuse this distance is 250 fis 
and 25 /is for wildtype and tail-deleted mutant, respec- 
tively. For smaller times, t < t c , k on plateaus at its initial 
value. For larger times, t > t c , it decays rapidly towards 
zero. The single molecule behavior is governed by the 
dimensionless number k = k on s 2 / AD , which is the ratio 
of timescales set by diffusion and rebinding. Diffusion 
does not interfere with rebinding as long as k > 1. Our 
theory therefore predicts that for wildtype with diffusion 
constant D = 10~ n cm 2 /s and capture radius s = 1 nm, 
k on > 4 x 10 3 Hz. For the tail-deleted mutant, mobility 
does interfere with rebinding and we must have k < 1. If 
we assume that in this case D is smaller by one order of 
magnitude, then k on < 4 x 10 4 Hz. Thus we can conclude 
that k on should be of the order of 10 4 Hz. 

Tether stabilization through multiple bonds. 

We now turn to the possibility that spatial proximity re- 
quired for rebinding is established by multiple contacts. 
Tethers above the shear threshold are modeled as clus- 
ters of N bonds, which in practice are expected to be 
distributed over at least two microvilli. At any time- 
point, each of the N bonds is either closed or open. The 
way force is shared between the closed bonds depends 
on the details of each tether realization. However, we 
expect that only those realizations will contribute signif- 
icantly to the long-lived tethers above the shear thresh- 
old in which different bonds share force more or less 
equally. This most likely corresponds to two microvilli 
being bound with similar latitude in regard to the direc- 
tion of shear flow. With this assumption, the force used 
in the single molecule dissociation rate has to be over- 
all force divided by the number of closed bonds. If one 
bond ruptures, force is redistributed among the remain- 
ing bonds. Open bonds can rebind with the rebinding 
rate k on . If rebinding occurs, force again is redistributed 
among the closed bonds. In general, in the absence of 
diffusion cluster lifetime T, but not the full cluster dis- 
sociation probability function can be calculated exactly 
|22|. We first discuss the case without loading or diffu- 
sion, thus focusing on the role of rebinding. As argued 
in the supplemental material, for small rebinding rate, 
k n < fco j cluster lifetime T scales logarithmically rather 



5 




0.2 0.4 0.6 

pause duration (s) 



FIG. 4: Computer simulations show that L-selectin mediated 
tethers above the shear threshold yield first order dissociation 
kinetics. Solid lines: two-bonded tether with F = and 
k on = 10 4 (right) and 0.5 x 10 4 Hz (left). Dashed lines: the 
same with F = 100 pN. Dotted lines: k on = 10 4 , F = and 
mobility parameter k = k on s 2 /AD — 1 (right) and 0.5 (left), 
respectively. Inset: L-selectin mediated tethers show Bell- like 
shear dependence even in the presence of L-selectin mobility. 
Solid lines from bottom to top: no mobility, k = 2.5, 1 and 
0.5. 



than linear with cluster lifetime N. This weak increase in 
T with N results because different bonds decay not one 
after the other, but on the same time scale. The exact 
treatment shows that for clusters of 2, 10, 100, 1000 and 
10, 000 bonds without rcbinding, lifetime is prolongued 
by 1.5, 2.9, 5.2, 7.5 and 9.8, respectively. In order to 
achieve 14-fold stabilization as observed experimentally 
at the shear threshold, one needs the astronomical num- 
ber of 6 x 10 5 bonds. In practice, for the case of di- 
lute ligand discussed here, only very few bonds are likely. 
Therefore even in the presence of multiple bonds, rebind- 
ing is essential to provide tether stabilization. 

In general, fast rebinding is much more efficient for 
tether stabilization than large cluster size. Our calcula- 
tions predict that in order to obtain 14-fold stabilization 
for the cases N — 2, 3 and 4, one needs k on = 6 x 10 3 , 
10 3 and 550 Hz, respectively. The value k on = 6 x 10 3 Hz 
obtained for the case N = 2 is surprisingly close to the 
estimate k on — 10 4 Hz obtained above via a completely 
different route, namely the competition of rebinding and 
diffusion for a single molecule. Therefore in the follow- 
ing we restrict ourselves to the simple case of two bonds 
being formed above the shear threshold (most probably 
by two microvilli) . In this case, cluster lifetime can be 
calculated to be |22| : 

T= J_( e -F/2F b+2e -F/F b + kon e -3F/2F b \ (1) 

2k V fa J 



A derivation of this result is given in the supplemental 
material. In Fig. 3 we use Eq. to plot the dissocia- 
tion rate for the two-bonded tether (identified with the 
inverse of cluster lifetime T) as a function of shear rate 
for different values of rebinding. The shear threshold at 
40 Hz corresponds to F = 0.36 Fb. It follows from Eq. Q 
that for this value of F, 14-fold stabilization in compar- 
ison with the force-free single bond lifetime is achieved 
for k on « 44 fa- For fa — 250 Hz, this corresponds to a 
rcbinding rate of k on = 1.1 x 10 4 Hz. Thus again we ar- 
rive at the same order of magnitude estimate, k on = 10 4 
Hz. Fig. 3 shows that with this value for k on , agreement 
between theory and experimental wild-type data above 
the shear threshold is surprisingly good. 

Relation to BIAcore. We now discuss how our esti- 
mate relates to BIAcore data for L-selectin [23j. In this 
experiment, L-selectin was free in solution and GlyCAM- 
1 immobilized on the sensor surface, which makes it a 
monovalent ligand. For the equilibrium dissociation con- 
stant, the authors found Kd = 105 /iM. This unusually 
low affinity results from a very large dissociation rate 
k r , which they estimated to be k r > 10 Hz. The re- 
sults presented in Fig. 1 seem to suggest that the real 
dissociation rate k r = 250 Hz. However, surface anchor- 
age of both counter-receptors often reduces bond lifetime 
by up to two orders of magnitude |24| . This has been 
demonstrated experimentally for several receptor-ligand 
systems and might result from the reduction in free en- 
thalpy of the anchored bond. Thus it might well be that 
the dissociation rate fa = 250 Hz found for surface an- 
chored bonds might be reduced down to k r = 10 Hz for 
free L-selectin binding to surface-bound ligand. Then 
the association rate kf = 10 5 1/Ms. For a capture ra- 
dius s = 1 nm and a three-dimensional diffusion constant 
D = 10~ 6 cm 2 /s, the diffusive forward rate in solution 
is k + = A-kDs = 8 x 10 8 1/Ms. Because k + > kf, the 
receptor-ligand binding in solution is reaction-limited, as 
it usually is. As explained above, k on can be identified 
with the rate with which an encounter complex tranforms 
into the final product 0, [^. Since bond formation 
is reaction-limited, k on = k/K + — 4 x 10 4 Hz, where 
K+ = 3/47r,s 3 is the dissociation constant of diffusion. 
Thus this estimate agrees well with the two other esti- 
mates derived above. 



Computer simulations 

In order to obtain effective dissociation rates in the 
presence of diffusion, one has to use computer simula- 
tions. For each parameter set of interest, we used Monte 
Carlo simulations to simulate 5,000 realizations accord- 
ing to the rates given above. More details are described 
in the supplemental material. In general, our simulations 
show that for strong rebinding, that is k on > fa, the effec- 
tive dissociation kinetics of small clusters is first order. In 
Fig. 4, this is demonstrated for the case N = 2. The plot 
shows the logarithm of the simulated number of tethers 
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lasting longer than time t for different parameter values 
of interest. All curves are linear, even in the presence of 
mobility, and the slopes can be identified with the disso- 
ciation rates. For example, k on — 10 4 Hz and F = 1 00 
pN yields the same effective first order dissociation rate 
as k on — 0.5 x f0 4 Hz and F = 0, thus rebinding can 
rescue the effect of force. Our simulations also show that 
cluster dissocation rate as a function of force fits well to 
the Bell equation for k on > ho- In particular, this holds 
true in the presence of L-selectin mobility, as shown in 
the inset of Fig. 4. For k on = 10 4 Hz and without mo- 
bility (vanishing diffusion constant), lifetime at the shear 
threshold is 12-fold increased compared with single bond 
dissociation. With our estimate for wildtype mobility 
(k = k on s 2 /AD — 2.5), 10-fold stabilization takes place. 
For tail-deleted mobility (k = 0.25), only 1.5-fold stabi- 
lization occurs. This effect is more dramatic than ob- 
served experimentally, where stabilization for wildtype 
and tail-deleted mutant are 14-fold and 7-fold, respec- 
tively. In practice, the mobility scenario is certainly more 
complicated and is expected to smooth out the threshold 
effect arising from our modeling. 



Discussion 

In this paper, we have presented biophysical modeling 
of L-selectin tether stabilization in shear flow based on 
recently published flow chamber data with high tempo- 
ral resolution . Our analysis suggests that the 14-fold 
stabilization observed at the shear threshold results from 
formation of multiple contacts and a single molecule re- 
binding rate of the order of k on — I0 4 Hz, which is re- 
markably faster than the force-independent dissociation 
rate fco = 250 Hz observed below the shear threshold. 
Using computer simulations, we showed that for such 
strong rebinding, the experimentally observed first order 
dissociation and Bell-like shear force dependence follow 
from the statistics of small clusters of bonds. Despite the 
good quantitative agreement achieved here between ex- 
perimental data and our model, it is important to state 
that it cannot be expected to predict all details of the ex- 
perimental results. In practice, the formation of bonds is 
a stochastic process and there will be a statistical mixture 
of differently sized and differently loaded clusters, involv- 
ing different microvilli and different scaffolds of L-selectin 
ligands. Cytoskeletal anchorage of the different ligand- 
occupied L-selectin molecules might also change in time 
in a complex way. Nevertheless, by focusing on the case 
of two bonds (possibly on two different microvilli) with 
shared loading and mobility-dependent rebinding, we ob- 
tained quantitative explanations for many conflicting ob- 
servations from flow chamber experiments and biomem- 
brane force probes, which have not been interpreted in a 
consistent way before. 

Several explanations have been proposed for the shear 
threshold effect before. Chang and Hammer suggested 
that faster transport leads to increased probability for re- 



ceptor ligand encounter |25j . Yet the new high resolution 
data from flow chamber experiments indicate that below 
the shear threshold, the issue is insufficient stabilization 
rather than insufficient ligand recognition [I5| . Chen and 
Springer suggested that increased shear helps to over- 
come a repulsive barrier, possibly resulting from negative 
charges on the mucin-like L-selectin ligands j4|. However, 
Dwir and coworkers showed that small oligopeptide lig- 
ands for L-selectin presented on non-mucin avidin scaf- 
folds exhibit the same shear dependence as their mucin 
counterparts |l5j . Evans and coworkers have argued that 
increased shear leads to cell flattening and bond forma- 
tion [l2| . However, Dwir and coworkers found that fix- 
ation of PSGL-1 presenting neutrophils does not change 
the properties of tethers formed on low density immobi- 
lized L-selectin, while they do destabilize PSGL-I teth- 
ers to immobilized P-selectin (Dwir and Alon, unpub- 
lished data). These data suggest that cell deformation as 
well as stretching and bending of microvilli do not play 
any significant role in L-selectin tether stabilization. Re- 
cently, the unusual molecular property of catch bonding 
has been suggested as explanation for the shear threshold 
[26l |27| . However, the data by Dwir and coworkers sug- 
gests that force-related processes do not account for the 
shear threshold of L-selectin mediated tethering [l5j . Our 
interpretation of the shear threshold as resulting from 
multiple bond formation is supported by experimental 
evidence that increased ligand density both rescues the 
diffusion defect and abolishes the shear threshold [l5| . 
The diffusion defect can also be rescued by anchoring 
of cell-free tail mutants of L-selectin to surfaces, allow- 
ing them to interact with leukocytes expressing L-selectin 
ligands |14| . 

On all ligands tested, the tail-truncated and more so 
the tail-deleted L-selectin mutants support considerably 
shorter tethers, consistent with a role for anchorage in 
these local stabilizaton events. One possible explana- 
tion is that cytoskeletal anchorage prevents uprooting of 
L-selectin from the cell. However, uprooting from the 
plasma membrane of neutrophils has been shown to take 
place on the timescale of seconds 28] . The tail-truncated 
L-selectin mutant has still two charged residues in the 
tail, which makes it impossible to extract it from the 
membrane in milliseconds. Receptor uprooting from the 
cytoskeleton only should lead to microvillus extension, 
which however is a slow process and has been shown 
to stabilize the longer-lived P-selectin mediated tethers 
rather than L-selectin mediated tethers |2j|. Here we 
postulated another possibility for cytoskeltal regulation, 
namely restriction of lateral mobility. It has been argued 
before for integrin-mediated adhesion that increased re- 
ceptor mobility due to unbinding from the cytoskeletal 
is used to upregulate cell adhesion [2(1 • Indeed in- 
creased receptor mobility is favorable for contact forma- 
tion, but here we show that it is unfavorable for contact 
maintenance, since it reduces the probability for rebind- 
ing. 

Our analysis suggests that the smallest functional teth- 
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ers are mediated by a least two L-selectin bonds, each 
on a different microvillus, working cooperatively as one 
small cluster. Our model does not explain from which 
configuration a broken bond rebinds, but it suggests that 
this configuration is neither collapsed (otherwise rebind- 
ing, which implies spatial proximity, was not possible) 
nor strongly occupied (otherwise diffusive escape was not 
possible). We can only speculate that complete rupture 
is a multi-stage process, and that the rebinding discussed 
here starts from some partially ruptured state. We also 
cannot exclude that the rebinding events described here 
involve different partners than the dissociated ones, be- 
cause both L-selectins and their carbohydrate ligands 
might be organized in a dimeric way. Moreover, cytoplas- 
mic anchorage might proceed in multiple steps, includ- 
ing some weak pre-ligand binding anchorage, which is 
strengthened by L-selectin occupancy with ligand. Cou- 
pling between ligand binding and cytoplasmic anchorage 
is well-known for integrins |30j and might also be at work 
with selectins. 

The mechanisms discussed in this paper could be effec- 
tive also with other vascular counterreceptors specialized 
to operate under shear flow. As argued here, the excep- 
tional capacity of L-selectin to promote functional adhe- 
sion in shear flow might not only result from fast dissoci- 
ation and high strength under loading, but more so from 
a fast rebinding rate. Indeed other vascular adhesion re- 
ceptors specialized to capture cells share on-rates similar 
to that of L-selectin [31|. Shear flow may also promote 
multi-contact formation for shear-promoted platelet teth- 
ering to von Willebrand factor [32] • It may also enhance 
formation of multivalent a40j and LFA-1 integrin tethers 
to their respective ligands [Ul^l- The importance of cy- 
toskeletal anchorage in local rebinding processes of these 
and related adhesion receptors has not been experimen- 
tally demonstrated to date. However, the lesson drawn 
here from the role of L-selectin anchorage in millisecond 
tether stabilization may apply to these receptors as well. 
Future studies will help confirm this hypothesis. They 
will also shed light on the specialized structural features 
acquired by these receptors and their ligands through 
evolution, allowing them to operate under the versatile 
conditions of vascular shear flow. 

Acknowlegdments: We thank Oren Dwir, Thorsten 
Erdmann, Evan Evans, Stefan Klumpp, Rudolf Merkel, 
Samuel Safran and Udo Seifcrt for helpful discussions. 
R.A. is the Incumbent of The Tauro Career Development 
Chair in Biomedical Research. U.S.S. is supported by the 
German Science Foundation through the Emmy Noether 
Program. 



Supplemental material 

Model. We consider a cluster with a constant number 
N of parallel bonds under constant force F. At any time 
t, i bonds are closed and N — i bonds are open (0 < i < 
N). Closed bonds are assumed to rupture according to 



the Bell equation |Toj : 

k of} = k Q e F ' F » . (2) 

For convenience, we introduce dimensionless variables: 
dimensionless time r = fcot, dimensionless dissociation 
rate fc ///fco and dimensionless force / = F/Fb- The i 
closed bonds are assumed to share force / equally, that 
is each closed bond is subject to the force f/i. Thus the 
dimensionless dissociation rate is e^ 1 . As long as the 
receptors are held in proximity to the ligands, rebinding 
of open bonds can occur. Therefore we assume that single 
open bonds rebind with the force independent association 
rate k on . The dimensionless rebinding rate is defined as 
7 = k on /k . 

The stochastic dynamics of the bond cluster can be 
described by a one-step Master equation [22| 

dpi 

— = n+ipi+i + g l -xv%~\ ~ [n + 9i\Pi (3) 
ar 

where Pi(r) is the probability that i closed bonds are 
present at time t. The reverse and forward rates between 
the different states i follow from the single molecule rates 
as 

n = ie"' , gi = 7 (JV - i) . (4) 

Once the completely dissociated state i = is reached, 
the cell will be carried away by shear flow and the cluster 
cannot regain stability. This corresponds to an absorbing 
boundary at i = 0, which can be implemented by setting 
r a = 9o = 0. Cluster lifetime T is identified with the 
mean time to reach the absorbing state i = 0. 

Lifetime of two bonded cluster. For a cluster with 
two bonds, N — 2, cluster lifetime T can be calculated 
exactly in the following way. At r = 0, the cluster 
starts with the initial condition i = 2. Next it moves 
to state i = 1 with probability I, after the mean time 
l/ r 2 ■ From there, it rebinds to state i = 2 with prob- 
ability wr, = gi/{r\ + gi), or dissociates with probabil- 
ity Wd = r\j(r\ + <7i). The mean time for this part is 
l/(ri +gi). Thus after two steps the system has reached 
state i = with probability wd or returned to state i = 2 
with probability wr, with wd + wr = 1- In detail, the 
probabilities and mean times for both processes are 

e f II 
w D = -r- — , t D = —jj^ + —r- — , (5) 

7 

wr = —7- , t R = t D . (6) 

e J + 7 

Different paths to dissociation only differ in the number 
of rebinding events j to state i = 2: 

wj = w D w R , tj =t D + jt R . (7) 

We first check normalization: 
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and then calculate cluster lifetime: 



T = ^ t 3 W 3 =t D + tRW D ^ 3 



3=0 

= t D + t R W D 



3=0 

WR t D 



(9) 
(10) 



(1 - WR) 2 1 - WR 

i( e -//a + 2e-' + 7 e" 3 - f / 2 ) . (11 



This formula is given in dimensional form as Eq. 1 in the 
main text. 

Cluster size versus rebinding. For arbitrary clus- 
ter size N, cluster lifetime T can be obtained from the 
adjoint Master equation [22L l35j. In the case of vanish- 
ing force, / = 0, the solution can also be found by using 
Laplace transforms [22|: 



T = 



(12) 



For 7 = 0, this equation reduces to 

N 



' * 1 



(13) 



i=l 



where are the harmonic numbers. An expansion for 
large N gives 



H 



N 



r + iniv 



i 

2N 



0( 



N 2 ' 



(14) 



Here T — 0.577 is the Euler constant. This formula is 
rather good already for small values of N. Eq. (|13J) is 
easy to understand: for 7 = 0, dissociation is simply a 
sequence of Poisson decays with mean times 1/rj = X/i. 
The overall mean time for dissociation is the sum of the 
mean times of the subprocesses. We conclude that for 
vanishing rebinding, T grows only weakly (logarithmi- 
cally) with N and very larg e cluster sizes are required 
to achieve long lifetimes [jj], |3(| . Therefore rebinding is 
essential to achieve stabilization for small cluster sizes. 

Effect of finite loading rate. Loading and dissoci- 
ation of single L-selectin bonds occur on the same time 
scale. As a cell is captured from shear flow and comes 
to a stop, force rises from cero and plateaus at a finite 
value. We model the initial rise as linear, with loading 
rate r. Therefore / = fXT until time To, followed by con- 
stant loading / = /o, where \x — r/koFb is dimensionless 
loading rate. Since fx = /o/to, there are only two inde- 
pendent parameters, tq and /o- The mean lifetime can 
be calculated in the usual way 0, 0] . We find 



T= — (e(-)-E{ ) + ^e-— 



(15) 



where E(x) is the exponential integral. For tq — > 0, 
we find the result for constant loading, T = 1/e-™. 
For tq — > 00, we find the result for linear loading, 
T = e 1/ ^E(ll 'n)/ 'n 17]. Eq. lO is used in the section on 
single bond loading and for the plot of the dash-dotted 
line in Fig. 3. 

Simulations. In the presence of diffusion with dif- 
fusion constant D, the single molecule association rate 
becomes a function of the time t which has passed since 
unbinding. In this paper, we use the approximation 



k n{t) — k on ^1 e I J 



(16) 



where s is the capture radius. Since anaytical solutions 
are intractable in this case, the Master equation Eq. © 
has to be solved numerically. The standard method to 
do so are Monte Carlo simulations. Unfortunately, the 
Gillespie algorithm for exact stochastic simulations [37j 
cannot be used in this case, because it does not track the 
identity of different bonds [33 ■ Therefore we simulate the 
Master equation by discretizing time r in small steps At. 
For each time step, random numbers are drawn in order 
to decide how the system evolves according to the rates 
defined for the different processes. In detail, in the time 
interval [r, r + At) each closed bond has the probability 
e l 7* At to rupture, and each open bond has the probabil- 
ity 7(1 - e~ fe/(7r) )Ar to rebind. Here k = k on s 2 /AD is 
the dimensionless ratio of the timescales set by diffusion 
and rebinding. Our model for stochastic cluster dynam- 
ics was implemented in the programming environment 
Matlab. A typical run simulates 5, 000 tethers (larger 
tether numbers give better statistics for the long time be- 
haviour, but similar results), each comprising N bonds. 
Results from different runs are bined into histograms for 
the number of tethers dissociating in the time interval 
[t, r + At]. In Fig. 4, we plot the natural logarithm 
of the fraction of tethers which last longer than dimen- 
sional time t as a function of t, as it is common for the 
analysis of experimental data. The slope of this curve is 
identified with the dissociation rate. Although this pro- 
cedure involves numerical integration of the probability 
distribution for dissociation, and therefore leads to loss 
of information, its smoothing effect is essential to obtain 
reliable estimates for the dissociation rate in the presence 
of noisy data. In the inset of Fig. 4, the dissociation rates 
obtained in this way are plotted as function of shear rate 
(that is force) and diffusion constant (which determines 
the dimensionless parameter k). 
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